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CN ' We discuss the implementation and properties of the quenched approximation in 

■ the calculation of the left-right, strong penguin contributions {i.e. Qq) to e'/e. The 



coefficient of the new chiral logarithm, discovered by Golterman and Pallante, which 
appears at leading order in quenched chiral perturbation theory is evaluated using 
both the method proposed by those authors and by an improved approach which is 
free of power divergent corrections. The result implies a large quenching artifact in 
Q I the contribution of to e'/e. This failure of the quenched approximation affects 



only the strong penguin operators and so does not affect the Qs contribution to e'/e 



^ . nor Re^O) Rej42 and thus, the AL = 1/2 rule at tree level in chiral perturbation 

H ; 

5r , theory. 
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I. INTRODUCTION 

There have been several recent apphcations of lattice QCD to the calculation of Re(e'/e), 
the direct CP violating parameters in K it n decavs. These include the attempts using 
domain wall fermions by the CP-PACS and RBC collaborations. A notable feature of 
both of these calculations is that their central values differ drastically from experiment. The 
experiments at CERN j^, 0] and Fermilab 0, 0] have yielded an experimental world average 
ofRe(eVe) = (1.8±0.4) x lO'^ fj. The RBC collaboration reported a value -4.0(2.3) x 10-^ 
and a similar negative central value was 

using staggered fermions and for ongoing work also using staggered fermions.) The stated 
errors in the RBC lattice calculations for e'/e were statistical only, with an estimate of the 
size of the systematic errors requiring much further study, of which this work is part. 

A number of serious approximations were made in the lattice calculations which intro- 
duced uncontrolled systematic errors. One of these was the quenched approximation, in 
which the fermion determinant is ignored in the generation of the gauge configurations. 
This is a truncation of the full theory that dramatically reduces the computer resources 
required, but is uncontrolled. Where quenched lattice results have been compared to ex- 
periment, in simple quantities such as masses of flavored mesons and decay constants, the 
agreement is at or better than ~ 15%. However, there is no apparent reason for this level 
of agreement to hold for all low-energy hadronic phenomena. Another approximation made 
in all existing lattice calculations of e'/ e was the use of leading order chiral perturbation 
theory (ChPT) to relate unphysical K —>■ it and K ^ amplitudes to the physical K tttt 
amplitudes, as first proposed by This is also likely to be a serious source of systematic 
error, although in this paper we focus on a particular ambiguity present in the quenched 
approximation. 



Since the original lattice calculations of |l[ Q , it was shown in 



12l | that, at leading 



order in quenched chiral perturbation theory there is a term logarithmic in the pion mass 
which contributes to the matrix elements of the strong penguin operators. This term is ab- 
sent in full QCD, and its contribution is proportional to an a priori unknown, new low energy 
constant (LEC). In terms of a representation of quenched QCD in which the fermion loops 
are cancelled by the addition of ghost fields to the lagrangian, this LEC can be associated 
with the presence of additional operators in the effective hamiltonian mediating the weak 
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decay of ^ vrvr which contain both quark and ghost fields. The presence of additional 
operators naturally calls for a re-examination of the way in which the quenched approx- 
imation is implemented in such matrix element calculations; without further guidance, a 
single strong penguin operator in the full theory can be represented by any arbitrary linear 
combination of its direct transcription into the quenched theory (a four-quark operator) and 
these additional (two-quark, two-ghost) operators. 

Since this ChPT result only became known after the original RBC analysis was completed 
P], that work used the leading order ChPT relevant for full QCD. As such, it is useful to 
discuss how the presence of these new operators might effect this result. We emphasize 
that this particular quenching difficulty is present only for the strong penguin operators 
and so does not affect the Qs contribution to e'/e. This quenching ambiguity also does 
not significantly affect Re Aq and Re A2 and since Re Aq receives only a negligibly small 
contribution from Qq Q], the RBC result for the AI = 1/2 rule remains unchanged. A 
recent paper penguin effects 

can arise in the AI =1/2 rule, and thus change the values of the underlying LEC's that one 
measures. However, the effects of Ref. [31 do not change the tree-level results of Bernard, 
et al 10], and as the original RBC analysis was performed at tree level, this does not affect 
the actual results for the A J =1/2 rule. An attempt to determine the systematic error due 
to quenching for this quantity may, however, be influenced by the results of Golterman and 
Pallante. 

Since the quenched approximation is uncontrolled, a rigorous matching of operators be- 
tween the quenched and full theories is not possible. However, we argue that the coefficients 
of these new operators can be determined by the same style of physical argument that is 
usually put forward to motivate the quenched approximation. This approach, which might 
be called "intermediate-energy matching" can be described as follows. Since quark loops 
play an important role in the renormalization group evolution of the weak amplitudes from 
the scale of the W, Z and top quark down to the kaon mass, the quenched approximation 
must be applied in a discriminating fashion, simulating the vacuum polarization of quark 
loops in the low-energy lattice QCD calculation with a weakened bare coupling while leaving 
the quark loops that appear in the perturbatively computed Wilson coefficients intact. 

Such a separation between short- and long-distance vacuum polarization effects can be 
made quite precise in the context of the effective weak Hamiltonian and gives a prescription 
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for carrying out the quenched calculation which is close to that adopted by RBC and CP- 
PACS. One simply requires at the intermediate energy scale at which the perturbatively 
determined weak amplitudes are matched to the low-energy four-quark operators in Hw 
that the full- and quenched-QCD amplitudes agree. Since there are no ghost quarks in 
the full theory, this requires that the ghost quark matrix elements of Hw evaluated in the 
quenched theory also vanish. This gives a physically motivated definition of the quenched 
approximation for the energy scale at which this matching is performed. Of course, as 
the quenched and full theories do not have the same low energy limit, physical results will 
depend on this matching scale; for the quenched approximation to be useful, the difference 
in these results as the matching scale is varied over the range of energy scales important 
for the quantities we are calculating must be numerically small. The extent to which this 
condition is obeyed depends on the quantity calculated; later we argue that for the strong 
penguin operators, the quenched approximation is particularly bad. 

Another possible guide we may take in transcribing the full QCD operators into 
the quenched theory is the extended chiral symmetry of this theory 

HQ. 

Recall 

that when the ghost quarks are added to normal QCD, the original chiral symmetry, 
SUl{3) X SUfi{3), which transforms only the normal quarks, expands into a larger, graded 
symmetry, SUl{3\3) x SUji{3\3) for the simplest case of three quarks u, d and s and three 
ghost quarks m, d and s. 

As analyzed by Golterman and Pallante, the original operators appearing in Hw-, trans- 
forming in specific representations of SUl{3) x SUr{3), take on new SUl{3\3) x SUr{3\3) 
quantum numbers. The new operators containing ghost quarks might then be chosen in a 
fashion to simplify the resulting representation of the extended SUl{3\3) x SUr{3\3) symme- 
try group while leaving the physical SUl{3) x SUr{3) behavior unchanged. However, since 
this extended SUl{3\3) x SUr{3\3) group is necessarily unphysical it is difficult to provide 
a convincing motivation for the mixture of ghost quarks chosen. 

Thus, we beheve that the quenched approximation can be applied to weak decays in a well- 
motivated way. One determines by "intermediate-energy matching" , the quenched effective 
Hamiltonian, H'^, and then examines matrix elements of the resulting operators analytically 
for possible quenched chiral logarithms and numerically to make quenched predictions from 
the theory. The size of the quenched chiral logarithms should be viewed as a measure of the 
errors in the quenched approximation. 
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From this perspective, the study of the quenched chiral logarithm discovered by Goiter- 
man and Pallante and the new low energy constant, a^'^ , which appears as its coefficient, is 



of central importance. In fact, these authors have proposed jl^ a method for obtaining a 
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directly from a lattice calculation. We show that there are difficulties using their method to 
obtain a^^ from lattice data using domain wall fermions, due to the presence of power diver- 
gences. However, motivated by their direct approach, we have found an improved method 
for obtaining a^^ that does not have these divergences. This is done by constructing an 
extension of the usual CPS symmetry arguments to the quenched case where both quark 
and "ghost" quark degrees of freedom are present. 

Using our proposed method, we obtain a value of a^^ which is consistent with the value 



lJ,lia. The value 



given by the large Nc approximation obtained by Golterman and Peris 
obtained from the method proposed by Golterman and Pallante also yielded a value roughly 
in agreement with the others, but with a larger systematic error. This large value of a^"^ 
has two important consequences: i) Large, non-analytic behavior in the quenched chiral 
limit which is absent in full QCD provides clear evidence of substantial systematic errors 
associated with the quenched approximation, ii) If these new quenched non-linearities are 
omitted from the functional forms used to extract the quenched LEC's, even the analysis 
within the quenched approximation is likely to be incorrect. 

In Section |n] we discuss the application of the quenched approximation to weak decay 
amplitudes and motivate the "intermediate-energy matching" approach described above. 
The chiral symmetry of partially quenched QCD and quenched chiral perturbation theory 
is reviewed in Sec. IIIII We follow the approach of Ref jlll and discuss the ambiguities of 
determirimg from this perspective. In Section HVl we discuss Golterman and Pallante's 



method 



12| for obtaining the low energy constant, cx^^ ■ We extend the usual CPS symmetry 



arguments for obtaining the form of the power divergent contributions to A J =1/2 matrix 
elements to the quenched case involving external ghost quark states. In Section El we use 
this extended symmetry, as well as numerical results, to show that the method of jl^ . 
when evaluated using domain wall fermions, suffers from ambiguities due to power divergent 
contributions. 

In Section IVII an alternative method for obtaining «f ^ is proposed which does not have 
a divergent contribution and this method is used in Section FVIII to obtain a numerical value 
for this constant. Finally, the implications of this result are discussed in the conclusion. 
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Section rVIIII Some useful formulae are given in Appendix A, our conventions are specified 
in detail in Appendix B, and a discussion of the one-loop ChPT calculation used in the 
numerical extrapolation to obtain a^^ is given in Appendix C. 

II. QUENCHED APPROXIMATION IN WEAK DECAY AMPLITUDES 

While the quenched approximation is widely used as a device to reduce the computational 
requirements of lattice QCD calculations, its implementation in the calculation of weak 
matrix elements deserves further discussion. The physical justification for this approximation 
is the hypothesis that the largest effect of the omitted quark loops is a modification of the 
running QCD coupling constant a<((/i) because of the omission of quark vacuum polarization. 
To the extent that this hypothesis holds true for a particular quantity of physical interest 
we can then compensate to a large extent for the effects of quenching on at a relevant 

physical scale /i by an appropriate weakening of the bare lattice coupling ao = Qq/^t^ where 
go is the gauge coupling which appears directly in the lattice QCD Lagrangian. An important 
limitation of this justification is apparent from the scale dependence of a<((/u). This scale 
dependence will be different between the full and quenched theories so the equality ^^(/i) = 
a'j.^{fi) will be approximately true only over a limited range of scales fi. 

Since a weak matrix element involves a large range of energy scales from the top quark 
mass down to Aqcd, a naive application of the quenched approximation to this entire energy 
range would be very inaccurate. However, because perturbation theory is used to describe 
most of this larger energy range, it is also not necessary to use the quenched approximation 
over the entire range. A clear view of this situation comes from considering some sample 
Feynman graphs representing the various ways that quark vacuum polarization loops can 
enter the gluonic penguin diagrams of interest. 

The possible effects of quark vacuum polarization are illustrated in Figs. QEl In the 
first of these. Fig. ^ the quark loop is entirely contained in the short distance part of the 
graph. The vacuum polarization loop will be dominated by momenta on the order of the 
top mass and can be accurately treated in QCD perturbation theory. Contributions with 
loop momentum on the order of Aqcd, where perturbation theory would not be accurate, 
will be suppressed by a factor of {Aqcd /^topY 

By contrast, the diagram shown in Fig. El contains a vacuum polarization loop that can 
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only contain low momentum. More precisely, once multiplicative charge and wave-function 
renormalization constants have been removed, any further contributions from large energies, 
for example fi ~ niw, will be suppressed by an additional factor of yU^/m^. Thus this graph 
represents a potentially non-perturbative piece which requires lattice techniques to evaluate. 
In fact, the above quenching hypothesis is simply the statement that the most important 
effect of removing this quark loop is a change in the charge renormalization, a change that 
can be completely compensated by a corresponding change in the bare lattice coupling g^. 

The diagram shown in Fig. El contains a vacuum polarization loop that can either enter a 
short- or long-distance part of the graph. These two possibilities are represented by the two 
dotted boxes shown in the figure. Because of the presence of the W propagator, the inner 
box is necessarily a "high-momentum" subgraph. However, important contributions can 
come from regions of momentum space in which the lines entering the vacuum polarization 
subgraph carry either small (~ Aqcd) ot large (~ mw) momentum. In the case that large 
momentum is involved, the outer dotted box surrounds what becomes a "high-momentum" 
subgraph with all internal lines carrying large momentum: a regime that can be evaluated 
perturbatively and one in which the quenched approximation should not be used. For 
the "low-momentum" case, the contribution will be non-perturbative and evaluation using 
lattice techniques, potentially using the quenched approximation, will be needed. Other 
momentum assignments for the lines in this outer box which include only a portion of the 
vacuum polarization subgraph will have more than four external quark lines, thereby having 
larger dimension and giving a contribution that is suppressed by a factor of (AQCD/''^iy)^- 

Finally, Fig.|3]shows a potentially more ambiguous case. Here, in addition to the possibil- 
ity that the vacuum polarization loop is contained entirely within a high- or low-momentum 
subgraph, (the outer and inner dotted boxes respectively) it may also be partially in both, as 
indicated by the middle dotted box. It is this intermediate subgraph which requires special 
discussion. For the case that all the momentum in the subgraph contained in the outer 
dotted box is large, the vacuum polarization loop can be treated perturbatively without 
recourse to the quenched approximation. For the case that only the momentum contained 
in the inner loop is large and those in the remainder of the graph are small, non-perturbative 
techniques will be needed and the quenched approximation may be used. Thus, in this case 
the vacuum polarization loop may be removed, or equivalently, a cancelling loop of ghost 
quarks included. 
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However, for the case of the intermediate dotted box which cuts the quark loop, the 
vacuum polarization piece is one-half within the high momentum part (the contents of this 
intermediate dotted box) but one half is to be evaluated in the low-momentum part. Here 
it is less obvious whether this loop is to be included perturbatively {i.e. incorporated in a 
Wilson coefficient) or to be cancelled by an added ghost quark contribution. In the language 
of the effective weak Hamiltonian H^, such a cancellation of this vacuum polarization graph 
would be achieved by including ghost operators in H^r. 

This variety of roles played by QCD vacuum polarization graphs might suggest that a 
precise definition of the quenched approximation for the evaluation of weak matrix elements 
would require new and elaborate development. However, as the above examples suggest, the 
standard field theoretic formulation of the "effective" low energy theory nicely deals with all 
of these questions, providing an unambiguous separation of the decay amplitudes into short- 
distance perturbative parts in which no quenched approximation is made and long-distance 
parts that must be evaluated non-perturbatively, possibly within a quenched approximation. 
The potential ambiguity associated with the inclusion of ghost operators is resolved by the 
matching conditions that are imposed to define the quenched effective theory. 

To be more concrete, consider the gluonic penguin portion of the effective low-energy 
weak Hamiltonian transforming in the (8, 1) representation of SU (3) l x SU (3)ij and written 
as a sum of four independent, dimension-six, four-quark operators: 

Hw = + C4Q4 + C5Q5 + CqQq (1) 

where {ci}i=s-6 are the four Wilson coefficients and {(5i}«=3-6 the four, conventional gluonic 
penguin operators: 

q=u,d,s 
q=u,d,s 
q=u,d,s 
q=u,d,s 

A sum over the repeated color indices a and b is understood. 

The above discussion of the contributions of the various vacuum polarization graphs to 
short- and long-distance physics mirrors closely the usual field-theoretic derivation of the 



-f)4g,7'^(l-7')g6 (2) 

- Y)d,q,r{l - ^')q, (3) 

-f)4g,7'^(l + 7')g, (4) 

-l')d,q,ra + l')qa. (5) 
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effective weak Hamiltonian in Eq. The usual "factorization" of weak amplitudes into 
short- and long-distance parts realized by provides exactly the needed separation of 
vacuum polarization effects into those that are correctly included using perturbation theory 
and those which are omitted in the quenched approximation — their effects being partially 
reproduced by a decrease in the bare coupling constant. 

Let us briefly recall this correspondence. The matrix elements of the local effective weak 
Hamiltonian given in Eq. ^ will accurately reproduce those of the complete theory if those 
matrix elements involve momenta which are small compared to the scale of the W meson 
and top quark masses — the scale at which non-trivial structure for the weak interactions 
becomes visible. 

Leading contributions in an expansion in l/m^^p and will come from regions of 

integration over internal momenta in which: a) All momenta in a particular subgraph are 
large, b) That subgraph contains any top quark and W boson internal lines, and c) That 
subgraph has itself the minimum number of external lines or, more precisely, represents 
the lowest possible mass dimension or largest possible degree of divergence. Under these 
circumstances, this subgraph can be treated as a structureless local composite operator made 
up of fields corresponding to the external lines of the subgraph. When evaluated, integration 
over this large-momentum region for the subgraph contributes to the coefficient Cj appearing 
in Hw 

In fact, the four coefficients {cj}i=3_6 can be simply defined as those required to make 
the complete and effective theories agree. Typically the Wilson operators {Qi}i=z-& will be 
defined by normalization conditions specified at an energy scale /i making both the operators 
Qi and coefficient functions Cj functions of this scale /i as well. Usually the scale /i is also 
the scale at which the Cj are determined by requiring specific Greens functions containing 
an Hw vertex to agree with those predicted by the complete theory. 

The application of this effective field theory formalism to the definition of the quenched 
approximation for the evaluation of weak matrix elements is now quite straight-forward. The 
vacuum polarization graphs which contribute to the short-distance/high-momentum part of 
the above analysis will necessarily contribute to the Wilson coefficients, can be evaluated in 
perturbation theory and need not involve the quenched approximation. Vacuum polarization 
effects which involve low momentum will enter the matrix elements of the Wilson operators 
Qi, will likely require non-perturbative techniques for evaluation and may be computed in 
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the quenched approximation. 

Now the equahty between matrix elements of the complete and the effective theories, 
imposed for four specific amplitudes at the scale /i will hold for a range of scales and a variety 
of matrix elements only to the extent that the quenched approximation is accurate and the 
matrix elements (e.g. {T[TT\Qi\K)) of the Qi in the unquenched theory agree with those 
evaluated in the quenched theory with the bare coupling weakened to compensate for the 
omitted quark-anti-quark screening. Note, in this definition of the quenched approximation 
we are replacing the matrix elements of the operators Qi, originally to be evaluated in full 
QCD, with matrix elements evaluated in a new theory: a theory in which ghost quarks have 
been introduced, following the procedures of Bernard and Golterman [16], to completely 
cancel the fermion determinant, and the gauge coupling at short distances has been decreased 
to compensate for this quark loop omission. 

Finally, let us examine the potential ambiguity associated with vacuum polarization loops, 
such as those in Fig. HI which contribute partially to the short- distance, perturbative Wilson 
coefficients and partially to the low-energy non-perturbative matrix elements. The "quench- 
ing" of the loop cut by the middle dotted box in Fig. |3] reduces to the question of whether we 
add to our four Wilson operators {Qi}i=z-& of Eqs. |21Elnew operators which contain ghost 
quarks. In fact, in the quenched theory there are four additional operators which have the 
same symmetry under SU{?>)l x SU{?>)r as these original four operators: 

Ql= Y._Sal'{l-l')da%l'{l-l'')qb (6) 

q=u,d,s 

Ql= J] ia7^(i-7'Hga7'^(i-7')% 

q=u,d,s 

Ql= ^.7^(l-7'Kl,7'^(l + 7')g6 

q=u,d,s 

Ql= 5^ ^.7^(i-7'Hg,7^(i + 7')ga 

q=u,d,s 

Including these operators in our quenched effective weak Hamiltonian will have the effect of 
introducing ghost quark loops which will (partially) cancel the problematic quark loop in 
Fig.H 

However, given the discussion above, the most consistent approach to the quenched ap- 
proximation is now clear. The effective, quenched weak Hamiltonian, should incorporate 



(7) 



(8) 



(9) 
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all eight possible operators: 

6 

1=3 

and the eight coefficients {cj}i=3_6 and {c*}j=3_6 should be chosen to make the matrix ele- 
ments of and those of the complete theory agree as closely as possible. This "matching" 
condition should be imposed at a sufficiently large scale fi that the required, perturbative 
evaluation of the complete theory is justified. On the other hand, given the inconsistent n 
dependence of quantities in the complete and quenched theories, the scale also should be 
chosen as close as possible to the low energy region in which the quenched matrix elements 
are to be evaluated to minimize these inconsistencies. The scale /i is thus an intermediate 
energy scale and this approach might be called intermediate energy matching. 

How are these four new coefficients {c*}j=3^4^5_6 to be evaluated and what values are 
expected? These new coefficients can be fixed by imposing conditions on the new ghost- 
quark amplitudes that appear in the quenched theory. They can be evaluated by comparing 
the complete theory, evaluated perturbatively, with the quenched theory evaluated either 
perturbatively or non-perturbatively. Since the complete theory does not contain ghost 
quarks, we should require that specific two-quark/two-ghost-quark amplitudes should vanish 
in the quenched theory. This would appropriately be done at off-shell, non-exceptional 
momentum at the scale fi and be imposed on color mixed and unmixed and left- and right- 
handed flavor-singlet combinations of ghost quarks. If carried out in perturbation theory, the 
absence of ghost quark coupling in the complete theory simply requires that all the c* vanish 
to leading order in a^. The choice c* = 0|j=3_4^5_6 is precisely the quenched approximation 
used by RBC and CP-PACS in their quenched kaon decay calculations, jjj, |^ . Because of 
the ghost quark couplings introduced by the self-contractions of Fig. these c* coefficients 
will be non-zero at order a^. While it would not be especially difficult to calculate the 
{c*}j=3^4^5_6 to order in perturbation theory or to evaluate them using the RI-MOM 
techniques employed for similar quantities by the RBC collaboration in Ref. P|, we expect 
that these effects will be quite small as were similar disconnected amplitudes that were 
evaluated in this earlier work. 

The majority of numerical simulations and the discussion above has been focused on 
the quenched approximation, in which all vacuum polarization loops are dropped from the 
calculation. However, there is additional insight and a useful framework for analysis to be 
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obtained in the case of "partial quenching" when a portion of the fermion determinant is 
included in the calculation. An important example involves the use of the wrong number 
of sea quarks, working with two rather than three flavors as in In fact with a positive 
Dirac operator (such as with domain wall fermions) and the proper algorithm one can 
simulate with a value of Nf varying continuously between and 3. The discussion for the 
quenched case applies quite directly to this partially quenched well. 

Again one attempts to approximate matrix elements computed in the full theory with 
those computed in the truncated theory without the full fermion determinant. The effects of 
the missing determinant are partially reproduced by adjusting the bare coupling to account 
for the missing polarization effects. Again the effective weak Hamiltonian to be used in a 
partially quenched theory should be chosen so that at an intermediate matching point, the 
full and the partially quenched theory agree. Clearly as Nf — > 3 this approximation will 
become increasingly accurate as the vacuum polarization effects of the two theories become 
identical. 

Even in the case Nf = 3 one can consider a partially quenched theory in which the 
valence quark masses do not match with the sea quark masses in the fermion determinant. 
For this case, the above discussion of quenching is simplified, since the differences between 
fermion masses appearing in quark loops can be neglected at the matching scale fi. Thus, 
the effective weak Hamiltonian, written in terms of the valence quark fields, will be the 
same as that in the full theory. Often setting mvai = ^sea, which can always be done, is less 
interesting than using mvai 7^ ^sea together with analytic results to explore the chiral limit. 

In this section we have proposed a definition of the quenched approximation which pro- 
vides a natural and self-consistent application of the standard quenched approximation to 
the case of hadronic weak decays. Up to corrections which are likely small, the calculations 
of both the RBC and CP-PACS collaborations use this definition in their evaluations of the 
gluonic penguin contributions, . However, since such an approximation is necessarily 
ad hoc and not systematic it must be used with suspicion. If unphysical effects, such as the 
quenched chiral logarithms discovered by Golterman and Pallante and evaluated numerically 
below, turn out to be large (as we will see they do), then the quenched approximation to 
these gluonic penguin amplitudes should be abandoned. 
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III. REVIEW OF STRONG PENGUINS IN QUENCHED CHPT 

Following Golterman and Pallante, we will now use chiral symmetry and chiral pertur- 
bation theory to study the effects of the quenched approximation on the gluonic penguin 
contributions to K meson decay. Chiral perturbation theory is an important tool which 
provides an approximation scheme in which two-pion decay amplitudes can be computed 
from vacuum and single-pion transition amplitudes. It can also provide an indication of the 
accuracy of the quenched approximation by identifying unphysical, quenched singularities 
which give rise to the infamous "quenched chiral logarithms" . 

In order to exploit chiral symmetry in the quenched approximation, we must adopt a 
field-theoretic_description of quenching. This can be done in two ways: the supersymmetric 
formulation )16] and the replica method Q]. As in the discussion above, we will adopt 
the original supersymmetric formulation and use the quenched chiral perturbation theory of 
Bernard and Golterman . In this method, the valence quarks are quenched by introduc- 
ing ghost quarks which have the same mass and quantum numbers as the valence quarks 
but opposite statistics. Therefore the ghost loops cancel the loops of the valence quarks, 
effectively setting the fermion determinant to a constant, which is precisely the quenched 
approximation. The chiral symmetry group of this action is SU{n\n)i x SU{n\n)R, a graded 
symmetry group, where n is the number of valence quark flavors. We adopt the definitions 
and notation of Bernard and Golterman needed for this supersymmetric approach (see 



also Ref. M). 

n 

As in Ref. j21|, one can also consider the partially quenched case, in which N sea quarks 
are introduced into the theory. In the notation of 21], one has n quarks, of which are sea 
quarks, so that there are n — N valence and, in addition, n — N ghost quarks. The valence 
quarks have arbitrary mass, while the sea quarks are chosen to be degenerate. In this case, 
the graded chiral symmetry group of the action is SU{n\n — N)l x SU{n\n — N)^. 



A. Notation and graded groups 

We will need to identify representations of the graded symmetry SU{n\n — N) group 
and will adopt the following notation. Since the fundamental representation of SU{N) and 
its complex conjugate are usually denoted, A^ and A^, we will adopt a similar description 
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of those representations of SU{n\n — N): {n\n — A^) and {n\n — A^). Thus, for example, 
a quark bihnear of the form Q{1 + 7^)Q', will belong to a representation easily identified 
as {{n\n — N)j^, {n\n — A^)_r). Here we use Q to represent a column vector whose first n 
components are the anti-commuting quark fields q and whose final n — N components are 
the ghost quark fields q. Likewise the quark bilinear Q'y^{l + 7^)Q' will belong to the 
product representation (l(n|n — N)l, {n\n — N)j^ x (n\n — N)ji) where we use the notation 
l{n\n — N) to identify the trivial representation of the group SU{n\n — A^). This represen- 
tation is easily constructed as the {2n — N) x {2n — A^) identity matrix /^'fe'""^'*' where the 
indices a and b transform as elements of the {n\n — N) and {n\n — N) representations of 
SU{n\n — N) respectively. This matrix has supertrace (defined below in Eq.lT^and Ref. [2^) 
str(/("l"-^)) = N. 

We will be interested in two irreducible representations which appear in the product 
{n\n — N) X [n\n — A^) above. The first is the trivial representation l{n\n — A^) already 
discussed. The second is the adjoint representation which we denote as adj(n|?T, — A^). This 
is the representation formed from the (2n — A^)^ — 1 generators ^ of the group SU{n\n — N). 
Here i identifies the element of adj(n|n — A^) while the indices a and b transform as elements 
of the {n\n — N) and {n\n — N) representations of SU{n\n — N) respectively. Since the 
matrices U defining SU{n\n — A^) have unit superdeterminant, sdet(?7)=l, the generators 
T* each have vanishing supertrace: str(T*) = 0. For the case A^ 7^ 0, the trivial and adjoint 
representations are entirely distinct. Operators belonging to l{n\n — N) cannot mix with 
those in adj(n|r;, — A^); the latter have a vanishing supertrace while the former do not. 

However, for the quenched case where there are no sea quarks, A^ = 0, their supertraces 
both vanish and an operator in the adj(n|n) representation can mix with a l{n\n) operator. 
In fact, in the quenched, A^ = case, the multiplication of the 2n fields q and g by a 
common phase factor has superdeterminant 1 making the unit matrix I^*^'"^ a valid generator 
of SU{n\n) which is no longer simple in this A^ = case. Thus, for the case A^ = the 
adjoint representation, adj(n|n), continues to have dimension (2n)^ — 1 but now includes the 
identity matrix /("I"). Of course there are still anti-hermitian matrices with non- vanishing 
supertrace. To include these we must extend the adjoint representation to a larger (2n)^- 
dimensional representation, denoted here adj(n|n) which includes the vectors in adj(n|n). 

Therefore, the product representation {n\n) x [n\n) must be decomposed in an unfamiliar 
way. In contrast with the A^ 7^ case, this representation is not reducible and cannot be 
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written as a direct sum of adj(n|n — A^) and l{n\n — N). The equation 

{n\n - N) X {n\n - N) = adj(ra|n - A^) + l{n\n - A^) (11) 

does not apply for the case N = 0. Instead this product forms the irreducible represen- 
tation adj(n|n). Contained within the (2n) ^-dimensional supervector space on which this 
representation acts is an invariant supervector subspace of dimension {2ny — 1 which forms 
the representation adj(?T,|n). Finally this supervector subspace of dimension (2n)^ — 1 itself 
contains a one-dimensional subspace which is again invariant under the original adj(n|n) 
representation matrices, transforming as the trivial l{n\n) representation: 

(n|n) X {n\n) = adj(?T.|n) D adj(n|n) D l(n|?T,) (12) 

The characteristics of this quenched, = case will be discussed further below. 



B. Quenched chiral symmetry of Qg 

Let us now examine the effect that the quenched approximation has on the chiral symme- 
try properties of the gluonic penguin operator Qq, the four-quark operator which is expected, 
along with Qs, to make the largest contribution to e'/e. Our discussion builds on that of 



Golterman and Pallante 



11 



12 1 . Recall that in the effective weak Hamiltonian for the full 



theory the Qq operator is given by Eq. repeated here for convenience: 

Qe = Sa7,{l - 7')4 Yl ^"^a + l')qa. (13) 

The right-hand factor in this operator is a sum over light flavors, q = u,d, s, so in the full 
theory this factor is a flavor singlet under the symmetry group SU As discussed in the 
previous section, in the (partially) quenched theory, vacuum polarization effects permit this 
operator to mix with a new operator which contains sea and ghost quarks and belongs to 
the singlet representation of SU{n\n — N)^- 

Ql = Sal,{l - i')d,J2Qbri^ + 7') 2a. (14) 

Q 

where the sum over Q contains all valence, sea and ghost quarks. An appropriate multiple 
of this operator can be subtracted from the original Qq operator, to create a new operator 
which transforms in the adjoint representation of SU{n\n — A^)r: 

Qf = Qe- ^Ql (15) 
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Here the factor of is easily chosen so that the adjoint operator Qg has a vanishing 
supertrace. Following Golterman and Pallante, we re-order Eq. El to express the original 
operator Qq in terms of Q'^-^ and : 

Qe = Qf + ^Ql (16) 

This is a useful equation because the adjoint representation, to which belongs, also 
includes the usual electro-weak penguin operator Qs so the matrix elements of these two 
operators are connected by a simple supersymmetry transformation. In addition, the NLO 
chiral perturbation theory for the relevant matrix elements of Qf has been worked out. 

As is evident from these equations, this useful decomposition does not work in the truly 
quenched case where iV = 0. In this case, Golterman and Pallante propose defining a 
non-singlet operator Q^^ through the equation: 

Qe = Qr' + iQf- (17) 

where we have added the extra Q to the superscript of the general operator Q| defined 
in Eq. El to emphasize that it is being defined for the quenched (A^ = 0) case. At first 
glance, Eq. El might offer the possibility of distinguishing two distinct contributions to 
Qq: due to Qq^^ and Qq^ , based upon their different chiral transformation properties 
under the graded symmetry. In turn this separation could - perhaps - be used to motivate 
an alternative definition of the quenched approximation, in which the non-singlet piece 
is dropped. Unfortunately, the decomposition in Eq. El is entirely arbitrary. The "non- 
singlet" operator, Q^^^ is not protected by the graded symmetry from mixing with the 
singlet operator. It is therefore ambiguous to separate out the "singlet piece" of Qq, any 
amount of Qq^ could be added to the definition of Qq^^ , allowing the coefficient of 1/2 in 
Eq. El to be replaced by any arbitrary number. In the representation theory language of the 
previous section, Q^^ is actually contained within the representation to which Q^^^ belongs: 
l{n\n) C adj(n|n). Thus, we cannot use Eq. Elto identify a possibly preferred "singlet" part 
of the original Qq operator. It should be noted that is protected by the graded symmetry 
from mixing with non-singlet operators; our conclusion due to the preceding arguments is 
simply that, since it is impossible to unambiguously define the singlet piece of Qq, any such 
approach is difficult to motivate physically. 

For concreteness, we have focused on the specific operator Qq. However, the transfor- 
mation properties of the other three gluonic operators are quite similar. The operator 
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transforms in an identical fashion as does Qq and these will of course mix when the energy 
scale at which they are defined is changed. The operators Q3 and Q4 are somewhat different 
since they transform only under SU{n\n — N)l. While this increases the number of repre- 
sentations that can appear, the absence of right-handed indices makes their quenched chiral 
perturbation theory less singular. 

We will now exploit these quenched chiral symmetry properties to study the matrix 
elements of these operators in the chiral limit. 



C. Review of quenched ChPT 



We now focus on the case of central interest in this work: quenched QCD, as above 
concentrating on the operator Qq. Since there are separate ChPT predictions for the matrix 
elements of the operators, and Qq we will analyze them both in this paper. To leading 
order in QChPT these operators can be represented by 



Qf 



a^i^str[A65^S9^S^] + 2 ag^^o str[A6(MS"f + SM^)] + h.c, 
i{afi'^^str[A6a^Ea'^S^] + 2 af/^^o str[A6(MS^ + SM^)] + h.c.} 
+af^str[A6SASt] +h.c. 



(18) 



(19) 



where (Ae)^^ = 5^3(5^2 while 



M 
A 
Bo 



rriu, rrid, rris, rriu, m^, m^jdiag, 

l)diag) 



1,-1,-1 

2 

m . 



m 



K+ 



m 



(20) 
(21) 
(22) 



rriu + md rriu + rris rud + rUs 

The matrix A is an element of the extended adjoint representation adj(n|n) and appears here 
because the corresponding operator in the underlying quenched theory, defined in Eq. El 
transforms in the same fashion. The meson field S is defined by 



exp 



/ 



with 




(23) 



(24) 
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The quantity / is the pseudoscalar decay constant in the chiral hmit (In the normahzation 
used the physical value of the pion decay constant is f^^ ~ 130 MeV) while the 3x3 matrices 
(p, (j) and cire constructed from Goldstone fields which create and destroy particles made 
from valence quarks and anti-quarks (bosons), ghost-anti-ghost quarks (bosons) and quarks 
and anti-ghost quarks (fermions) respectively. 

We note that the ChPT representation for the operator Qq given in Eq. differs from 
the implications of the original formula of Golterman and Pallante, Eq. 3.5 of Ref. In 
particular, the quenched non-singlet operator Q^^^ will be represented in chiral perturbation 
theory by both singlet and non-singlet operators. As a result, the low energy constants a)^{ 
and a)^2 which multiply the two quenched singlet operators which appear in Qq need not 
agree with the coefficients a^i and which appear in the singlet operator Q^^ . 

In Eqs. El and ^1 and those below, we combine the 2n x 2n matrix S together with the 
similar matrices A and M, which appear in the quark-level theory, to form the most general 
set of operators that are invariant under the complete graded symmetry S{n\n)L x SU{n\n)R 
of the quenched theory, to a given order in the Goldstone particle masses and momenta. As 



matrices transforming in the product representation (n|n) x {n\n) they can be written as: 




where the sub-matrices have the same dimension as the sub-matrices of $, above. The 
invariance under the graded symmetry of the quenched theory requires the presence of 
supertraces in the operators of Eqs. UHl and [121 defined by 

str(f/) =tr(A) -tr(D). (26) 

To tree-level in ChPT, a pseudoscalar meson mass is given by 



m. 



,2 



Bo{mi + mj), (27) 



where and m,- are the masses of the two quarks that form the meson. We define m^^ to 



be the tree-level meson mass of two valence strange quarks, as in 



2l| 



mi, = 2mj, - mi. (28) 

The new non-singlet operator of Eq. ^lis nominally of 0{p^) in ChPT, but its tree-level 
contributions to physical matrix elements vanish because A is proportional to the unit matrix 
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in the valence sector 12 1. At O(p^) this is not true and the one- loop insertions of Qq^^ 
make a contribution of the same order as the tree-level insertions of Q^^ ■ Of course, there 
are additional non-singlet operators that can be constructed from the matrix A which enter 
at O(p^) whose presence is needed to compensate for the scale dependence of the one-loop 
insertions of Qq^^ . This introduces three more LEC's into the amplitudes we must consider 
in this paper. The effective Lagrangian to this order is [l^ : 

^'Ns°' = j:^r'or', (29) 

i 

with 

C»P = strfAeL^StASL'^], 

03^^ = str[A6{StAS,L2}], (30) 
Of^ = str[A6{StAS,S}], 
and S = 2Bq{M'^I1 + S^M), = zS^S^S. In the following we work with operators renor- 



malized in the MS scheme, absorbing the divergence into the coefficients of the effective 
Lagrangian, which have the form 



NS ^ r,NS ^ 1 



+ ^(7£; - ln47r) 



2a^'Vi- (31) 



2 

The finite coefficients, the renormalized low energy constants of the theory, 

while the factors rji are chosen to cancel the divergences of the one-loop insertions of the 
tree-level operator. In the quenched theory we find r^i = 0, 773 = 3 and 774 = —3. The scale 
dependence of the LEC's is given to one-loop order by 

Ci (/^2) = c/ (/ii) + — — -In— , (32) 

[ATrfy /i2 

where /ii and /i2 are two different values of the chiral scale. The result for physical amplitudes 
should be independent of the scale and Eq. 1221 is obtained by requiring that all one-loop 
amplitudes for the Qq operator be scale independent when the O(p^) NLO LEC's are included 
in the calculation. 

The relations between our conventions for the LEC's and those of Golterman and Pallante 
11,^^ are given below. The constant a^^ has been chosen to have a normalization that 
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agrees with that of in Ref. 



0: 



a 



NS 
q,GP 



(33) 



In addition, we are working with the notation of Ref. [2^ for the NLO LEC's. Although that 
work dealt with the electro-weak penguins, both transform in the (adj(n|n — A^)^, adj(n|n — 
N)Ft) representation of the partially quenched graded symmetry in the case 7^ and we 
therefore keep that notation. The relationship between the two is 



= (4vr)^2cf ^ 
P^,' = (4vr)^2cr, 

P^,' = (4vr)^2cr. (34) 



D. Qq amplitudes in quenched ChPT 

We now review the ingredients necessary to obtain the contribution of both the operators 
Qq^ and Qq to the K nn amplitude to leading order, 0{p'^), in quenched ChPT. The 
chiral behavior of Qq is the same as that of in the full theory: 

4m«^ 



{n^n-\Qr\K') = ^{m^ - (35) 



where, as in the full QCD case, the needed LEG a^^ can be extracted from the — > and 
K ^ 7C matrix elements: 

mr\K') = ^-^{ml-ml), (36) 

{tt^\Q^^\K^) = -^a'^frriKm^ - ^a'^^TTiK. (37) 
The K TTTT matrix element of Qq includes the one-loop contributions of a^^ and the 
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OijP') LEC's, cf^. In this case, — vrvr is given by 



(vr+TT-IQel^") 



2 



(8,1) 



12 (m^ 



In 4 



ml- 9 
2—^ + 2m^ 



In 



m 



is: 



m: 



+ 



m 



4 ^ ^ 



m'' mt, 

TV 71 

2 / i7/_2 _2 _2 



+ lOm^ - 4m^ In 



ml 

2 /l™2 ^ 



+2m;^ F{m^, m^, -m^) - 2m9{mj^ - 4m: 




4mE 



m 



2m| 



(2F(m2 , m^, -ml) + F(m^, m^g, -m^)) 



+6m 



m 



ml 



(38) 



where the function F is given in Appendix A. Note the scale dependence in the logarithmic 
term proportional to a^^ ■ This scale dependence cancels that of c 
Eq. ins The Qq amplitudes for — > and K ^ ti are 



^•^5 as can be seen from 



(o|g6|i^°) 



(7r+|g6|i^+) 



4z 

7 



+ 2c: 



*g2 



TVS 
4 



m 



K 



m. 



TVS 



167r2 



m^ln^-2m^ln4+m^3ln^ 

/x^ /i^ /i"' 



-SCm^ - ml\ 



P 



1 (8,1) 



(8,1) 
2«.2 



2c: 



2c 



NS 



(39) 
(40) 



where in the expression for i^' — > vr we have set the quark masses to be degenerate (m^ 



m: 



m 



a 



q^'^"* and a^2^'' entering the equations for (Eqs. 



as in the numerical simulation |2]. Note, as discussed above, the singlet LEC's 

andlini) need not be the same as 



the corresponding LEC's, a^^ and a^'^ appearing in the corresponding expressions for Q 
(Eqs.EniinniandlSIl). 

Notice that the same linear combinations of LEC's appear in the above three equations. 



QS 



(8,1) 
•^91 



— 2c^^ and a^g'^V^ + 2c^'^. If the value of a^^ is small, then the procedure 



for obtaining K 



TTTT in which one neglects a^^ reduces to that of the full theory. In this 
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case, the previously mentioned combinations of LEC's replace a^i'^^ and 0:2^'^'', respectively. 
The strategy that was employed in Refs. Q and Q implicitly made this approximation. 



It has recently been suggested, however, that the va^ 
other LEC's in the amplitudes and cannot be neglected 
to obtain the result (in our conventions). 



ue of a^^ is large compared to the 



3 



A large Nc expansion was used 



«f = -\f'Bl (41) 

When quenched lattice values for these constants are substituted into this equation, one finds 
that the resulting value of a^^ is so large that it cannot be ignored. The non-linearities 
implied by the presence of a^^ in Eq. EHlneed to be included when this expression is used to 
extract the LEC's from the K —>■ lattice data and the explicit form of the a^'^-dependent 
term in Eq. EH needs to be taken into account in a quenched prediction for K tttt. Of 
course, with such large quenching artifacts, such a quenched prediction will be of limited 
value. 

It is therefore crucial to have a means for determining a^^ directly from a lattice calcula- 
tion. Golterman and Pallante have provided a method for doing just that. However, we have 
found that their method suffers from ambiguities which are power divergent in the lattice 
spacing when employed in the case of domain wall fermions, making accurate extraction of 
^Ns fj-QT^ ^i^Q lattice data rather difficult. In particular, with a finite Lg (the separation 
between the two physical, four-dimensional walls in the fifth dimension), this power diver- 
gent contribution, although suppressed by a factor of order of the residual chiral symmetry 
breaking, 0{m,.es), could be large. We discuss this difficulty in the following sections and 
present an alternative method for obtaining a^'^ which avoids this problem. 

IV. LATTICE DETERMINATION OF FROM K ^0 

In this section we review the method proposed by Golterman and Pallante to obtain a^^ 

n 

on the lattice [12], together with the (CPS) symmetry arguments needed to understand the 
form of the power divergences in the matrix elements of four-quark operators. As will be 
explained in the following, Golterman and Pallante introduced a new operator, Q^^'^, and 
extracted a^^ from matrix elements of this operator which included ghost particles in the 
external states. By extending the standard CPS symmetry arguments to include quark-ghost 
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transformations, we show that this method suffers from contamination from terms which 
are power divergent in the lattice spacing. Nevertheless, we present numerical results from 
this approach. While these results are in rough agreement with the large Nc estimates jl^ . 
due to the ambiguities associated with the power divergence the results must be considered 
inconclusive; in Section VI we suggest an alternative approach which does not suffer from 
this problem. 



A. Review of Golterman and Pallante's method 

In Eq. EHlit is difficult to numerically disentangle the logarithmic a^^ term from the linear 
term and possible higher order effects. It was for this reason that a new matrix element was 
suggested in Ref. 1^ to which a^^ contributes at 0{p^), so that it may be obtained more 
readily in the chiral limit. This can be done if one considers a matrix element where the 
ghost quarks can appear on external lines. In order to obtain such a matrix element, one 
must perform an SU{3\3)l flavor rotation of the operator Q^^^ into 



\ 1 q 
= -^QA67^(1-7')Q 2^7^(1 + 7')Q, (42) 

where now the d is a ghost quark field. The matrix \q is given by = Si^Sj^, a quenched 

chiral transform of the matrix Xq defined earlier. To leading order in ChPT, this operator is 



qQ^^ = «f^str[A6SAS^] +h.c. (43) 

Note that some care must be taken in order to maintain consistency in the sign conventions 
between the chiral and quark level operators. Our conventions are presented in detail in 
Appendix B. Since the above operator is in the same irreducible representation as Qq^^ , it 
is parameterized by the same low energy constants. Considering the matrix element K ^ 0, 
we have, to leading order 

mr'\K') = y«f . (44) 
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Although this method isolates at leading order, the NLO contribution has a power 
divergent coefficient, making the numerical extraction problematic. We discuss the way in 
which the mixing of four-quark operators with power divergent lower dimensional operators 
can be constrained by the symmetries of the theory in the following subsection. 



B. Power divergences and CPS symmetry 

In general, the A/ = 1/2 matrix elements of four-quark operators have a power divergent 
part, due to mixing with lower dimensional operators. This power divergence will involve 
the quark bilinears, 'sd and 575^, times a momentum independent coefficient One can 
define the following quark bilinear operator jlo| . 

e(3.3) = s(l-75)ci, (45) 
which is equal to a*^^'^)Tr(A6S) to lowest order in ChPT, where in our conventions, a'^^'^^ = 

^ rn 

We briefly review the use of CPS symmetry 23] to determine the form of the power 
divergences that enter our computations due to mixing with the lower dimensional operator 
of Eq. 1321 Here C and P are the usual charge conjugation and parity inversion symmetries, 
while S is the symmetry under exchange of the s and d quarks, which is exact when the 
quark masses are equal. 

The parity even part of the above operator, srf, has a CPS of +1, while the parity odd 
part, S75(i, has a CPS of —1. The operator, Qq^^ , has a CPS of +1 and the matrix element 
(OIQq'^^IK) is a parity odd transition. Therefore, the power divergence of this matrix 
element must be proportional to the matrix element of a lower dimensional operator which 
is also parity odd and has CPS=+1. We see that the only operator with these symmetries 
that can be constructed is ^756? multiplied by mg — m^. Thus, the power divergent part of 
{0\Qq'^^\K) is proportional to — m^. 

For the parity even transition, {7i\Qq^^\K), the mixing must be with a parity even lower 
dimensional operator with CPS= +1. In this case the only such operators are 'sd and 
{ms + iTLfij'sd. The former is ruled out by chiral symmetry. As domain wall fermions at finite 
Lg break chiral symmetry such a mixing can occur, albeit greatly suppressed due to the 
domain wall fermion mechanism. Thus, the dominant part of the power divergence of the 



25 



K ^ TT matrix element is proportional to nis + 771^. Since this is a statement at the operator 
level, it is true to all orders in the chiral expansion, an important result when dealing with 
the numerical data of a lattice calculation. 

We now wish to consider matrix elements of the operator Qq^^ , since this will give us 
a^'^ at 0{p^) in the chiral expansion. In order to proceed, we construct a new graded 
symmetry, CPS, where now S not only exchanges s and d quarks, but also exchanges all 
valence quarks with their corresponding ghosts, i. e. under S the six "quark" fields transform 
as 

(m, d, s, u, d, s) — > (tt, s, d, u, s, d). (46) 

In order to implement this CPS symmetry, we need to understand how C and P act on 
ghost quark states. 

Under C, the valence quarks have the following familiar transformation properties 

q ^ Clf, (47) 
q -q^C-\ (48) 

Thus, the quark bilinears, Id and S7^(i transform as, 

sd ds, (49) 
s'lf.d -d'y^s, (50) 

where we have used the identities, = and C^^'j^C = — 7J and the fact that the 
quark fields anti-commute. 

In the case of ghost quarks we can define C analogously to the case of valence quarks, 
such that 

q^cf. (51) 

However, the anti-ghost does not transform independently from the ghost, because it obeys 
bose statistics. Taking the transpose conjugate of both sides yields 



rC-\ (52) 



which differs in sign from Eq. Using the transformation rules in Eqs. |^andE21and the 
fact that ghost quark fields commute, we see that 

s'jfj.d -d'Ji.s, (53) 
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where the sign change demonstrates that this transformation for ghost quarks does indeed 
act as a charge conjugation. For bihnears that are made of a valence quark and an antighost 
quark the fact that the fields commute will cause a sign change compared to Eqs. EHl and 
IHUl For example, 

s-ff,d^ d'j^s. (54) 

Under P, the quark fields transform as 

q ^ Pq, (55) 
q qP, (56) 

where P^P = 1. It is simple to recognize that the same transformation under parity is also 
valid for the ghost fields. 

We are now ready to examine the CPS transformation properties of the relevant operators. 
The bilinear, sd, has a CPS of —1, while 5756? has a CPS of +1. The four quark operator, 
Q^^"^, has a CPS of +1. Thus, the matrix element, {0\Q'^^^\K) , which is a parity odd 
transition, must have a power divergence given by the matrix element of a lower dimensional 
operator which is both parity odd and has CPS=+1. We see that the bilinear operator s'j^d, 
which already has the same CPS as Qq^^ , can only have the coefficient + m^. This is 
opposite to the situation for the K amplitude, where the CPS of Qq'^^ and S'j^d 
was opposite, implying a power divergence proportional to — m^. The reason for the 
difference is the exchange of valence and ghost quarks in the right-handed part of Eq. 021 
under S, yielding an extra relative minus sign. 

Note, this same conclusion could also be reached by starting with the dimension-3 oper- 
ator which represents the quadratic divergence present in the usual operator Qq: 

q\M{l + 7^)g + qM^X&{l - 7^)g = (m^ + m,)srf + (m^ - m,)s7^d. (57) 

Note, the form of the left-hand-side of Eq.EZlis determined by the chiral symmetry properties 
of the matrix Ae (an element of (8, 1)) and the quark mass matrix M (an element of (3, 3)). 
The relative sign of the two terms on the left-hand-side of this equation is determined by 
CPS symmetry. The corresponding dimension-3 operator that mixes with Qq^''^ must have 
the same form except for the addition of the matrix A and a quenched chiral rotation of Ag 
to Ae to agree with the quenched chiral symmetry properties of Qq'^^: 

QXeMA{l + -f^)Q+QAM%{l - 7^)Q = (m, - md)sd- (m, + md)s-f^d. (58) 
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This argument demonstrates that K has a power divergence proportional to ms+rrid. 
When the fermion discretization does not preserve exact chiral symmetry (such as domain 
wall fermions with finite Ls), there will also be mixing with the operator 's'y^d. For domain 
wall fermions this is suppressed by a power of rrires- Such power divergences lead to large 
uncertainties when one tries to obtain the chiral limit of a matrix element without exact 
chiral symmetry on the lattice. However, these symmetry arguments suggest a way around 
this difficulty and we provide an alternative method for obtaining a^^ in Section IVll which 
makes use of the following observation: the matrix element {tt\Q'^^^\K) is parity even 
and since the parity even bilinear, sd, has CPS = —1, the divergence of i^' ^ tt must 
be proportional to — rrid. For degenerate quark masses, the power divergence vanishes 
completely. Again, this is a statement at the operator level and holds to all orders in the 
chiral expansion. 



C. Lattice contractions for Golterman and Pallante's method 



Continuing the derivation of Ref. 12], we carry out the Wick contractions for 



{^\Qpenguin\K^)^ where we specialize to the case, Q'^Zguin = Q^^^ , 

{OIQ^'IK') = Yl [lMiV,A,q)-lMiA,V,q)] (59) 

qS{u,d,s} 

+^ [il'^iV, A, d) - tl'^M, V, d) + tI'^{V, A, s) - iI'j^iA, V, s)] , 

with 

hi{j,k,q) = TTc{TTs[TjSd{xop,xo)Ss{xop,xoy-f5] 

xTi s[T k Sg{x op, X op)]}, (60) 

lMU,k,q') = TTs{TTc[TjSd{Xop,Xo)Ss{Xop,Xo)''-f5] 

xTi c[TkSg'{xop,Xop)]}. (61) 

In these contractions, Fy = 7^ and = 'jfj^'y^ and the traces are over spin or color. The 
quantity Sq{x,y) is the Dirac propagator connecting positions x and y for a quark of type 
q = [u, d, s). The position xq locates the source for the meson while Xop is the position of 
the operator Q^^^ . The derivation of Eq. EHl makes use of the fact that ghost and valence 
propagators are equal flavor by flavor (when properly ordered; see Appendix B). As discussed 
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in Ref. [1^, this is important because it allows us to obtain a^^ from contractions that have 
already been computed for {0\Qq'"^\K^) . 

Golterman and Pallante make the additional observation that the linear combination 
of contractions for {0\Qq^^\K^) is the same as Eq. EHl but with opposite sign for all but 
the I'j^^{j,k,d) terms. However, a careful treatment of the conventions for external ghost 
states shows that the K ^ matrix element has the opposite sign compared to that given 
in Golterman and Pallante [14], so that the contractions for {0\Qq^^\K^) are the same as 
Eq. but with opposite sign for only the I'mU, k, d) terms. Thus, since K ^ Q does not 
contribute at 0{jP), subtracting K ^ Q from i^T — > yields 

(0|Qr^|ir°) - (OlQr^li^-O) = i[WA,d)-I',,{A,V,d)] 

= jc^r + 0{p-). (62) 

Again, we point out that the OijP') terms are multiplied by a quadratic divergence in the 
lattice spacing. 

Using the leading order ChPT result, 

(0|G(3'3)|irO) = ^a(3,3) ^ _(o|s75rf|irO), (63) 
we obtain the following ratio: 

(64) 

Since we are dividing by precisely the matrix element that, according to the operator level 
discussion of the previous section, multiplies the power divergence, Eq.lMlexplicitlv gives the 
form of the power divergence, with the neglected higher order terms in the chiral expansion 
being free of these divergences. Again, a^^ can be obtained in the chiral limit. The power 
divergence of Eq. |S3]is proportional to m^, since it is the sum of two terms, one proportional 
to + 171(1 and the other one proportional to — frid, with the term cancelling between 
them. The only divergent term in this ratio depends linearly on rrid as shown in Eq. IMl 
However, after this term proportional to rrid has been removed by extrapolation to rrid — »• 
there remains an 0{mres/0''^) contribution for domain wall fermions, due to the power 
divergence and the lack of exact chiral symmetry for finite Q]. For our simulations, this 
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effect could be as large as the constant we are trying to extract, as we show in the next 
subsection. 

V. NUMERICAL RESULTS FOR a^^ FROM K 

We now examine the ratio of matrix elements of Eq. |M1 which can be obtained from the 
contractions given in Eq. E^and the value for {Q\'s'~^^d\K) . We briefly review the details of 
the ensemble used in Ref. to generate the needed contractions. 

The quenched configurations were generated with the Wilson gauge action at a coupling 
oi (3 = 6.0 with lattice volume 16^ x 32. The ensemble is comprised of 400 configurations, 
separated by 10,000 sweeps, with each sweep consisting of a simple two-subgroup heat-bath 
update of each link (Cabibbo-Marinari with Kennedy-Pendleton accept-reject step). The 
lattice cut-off was = 1.922(40) GeV set by the p mass. The domain wall fermion fifth 
dimension was Ls = 16 sites with a domain wall height M5 = 1.8. The resulting residual 
quark mass was 0.00124(5) in lattice units. For comparison, the value of corresponding 
to a pseudo-scalar state made of degenerate quarks with mass equal to the physical kaon at 
f3 = 6.0 is approximately 0.0185. The light quark masses in units of the lattice spacing were 
taken to be nif = 0.01, 0.02, 0.03, 0.04 and 0.05. 

The data for the ratio in Eq. IHU are given in Table HI and are plotted in Fig. EJ as a 
function of m^. We see from Eq. IHU that the intercept of this graph allows one to obtain 
a^'^ , modulo the effects of the residual chiral symmetry breaking. The two lines in Figure IHl 
represent linear fits to the data, the top one plotted as a function of m^, while the bottom is 
the same data plotted against md + rrires- While in the presence of residual chiral symmetry 
breaking effects neither of these two methods is known to be correct, comparing these two 
approaches allows us to estimate the order of magnitude of the 0(mres)/a^ ambiguity. As 
can be seen the difference is of the order of 30%, suggesting a potentially large uncontrolled 
error due to the power divergence. 

The result of an uncorrelated fit to the form of f]o + Vi{''T^d + ^res) "was: r^o = —5.94(14) x 
10-3^ r]i = -2.0781(33), with a xVdof of 1.67(18). The value of a^^ obtained from this fit 
was -1.20(3) X 10"^ where the error is statistical only. If one fits the data in Table H] to a 
linear function of (instead of + rrires), then the a^^ obtained is —1.72(3) x 10^^. For 
comparison, the large Nc approximation, Eq. |^ gives a value of a^^ = —1.63(48) x 10~^ 
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where, again, the error is statistical only. Here we have used the quenched lattice values 



2J|, obtained on lattices with the 



for / and Bq reported in the previous RBC works, 0, 
same size and coupling as those used here to determine ■ In Ref. j2J| it was found 
that / = 0.0713(53) and Bq = 1.59(3) in the same lattice units, for an earlier set of 85 
configurations. We see that there is rough agreement between the large Nc value and the 
range of values given directly by the lattice K ^ matrix element in the chiral limit. Since 
an a^^ of this size would cause significant changes in the K -kit amplitude for Qq, as 
mentioned previously, it is crucial to remove the large systematic error due to the residual 
power divergence. 

VI. LATTICE DETERMINATION OF ^ FROM K 

As we proved in Sec. IIVB| the matrix element, (vrlQ^^'^l-ft'), does not have any power 
divergences when the quark masses are degenerate. To NLO in ChPT, it is given by 

4 



{n^\Qr'\K^) = -i-X' + Oip')- (65) 



Again, we have an amplitude where a^^ can be obtained in the chiral limit, this time 



Q 



without any power divergent ambiguities. As discussed in the next section on numerical fits, 
the NLO logarithmic term vanishes for this matrix element in the degenerate mass case (See 
Appendix C for details of this calculation). The combination of contractions needed for the 
new amplitude (for degenerate quark masses) is 

+ [L'MiV,q)-Li,iA,q)], (66) 

q€:{u,d,s) 

with 

Llj{j) = Tr,{Trc[rjS'd(xop,Xi)S'„(xop,Xi)^75] 

xTTc[TjSu{Xop,Xo)Ss{Xop,Xo)''-f5]}, (67) 

LmU^q) = Trc{Tr4rj5'rf(a;op,a;i)75S'„(xi,Xo)5's(a;op,Xo)S5] 

xTi s[T j S g{x op, X op)]}. (68) 

Here we use the same notation as in Eqs. IHOlandl^and again exploit the fact that the ghost 



and valence propagators are equal flavor by flavor when properly ordered 
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VII. NUMERICAL RESULTS FOR a^^ FROM K 

The results for the {tx\Q'^^^\K) matrix element are given in Table. ITTl and plotted in Fig.[7| 
The chiral limit was obtained by a simple linear extrapolation of the form do + di{mf + mres)-i 
with do = 9.80(76) x lO"^, di = 0.624(17) and a xV^of = 0.10(11). This yields a value 
of a^^ = -1.24(10) X 10~^ In this result we neglect the small error in It turns out 
that the linear fit is exact to NLO in ChPT, as a direct one-loop calculation shows that in 
this degenerate mass case the logarithmic term vanishes. This calculation is discussed in 
Appendix C. 

A comparison between the two methods for obtaining a^^ , from K ^ n and from K 0, 
is shown in Figure where p/{2a^^'^^){n+\Q'^^^\K+) is plotted over the K results. 
With this normalization the K —>-tc matrix element given in EqlHSlhas the same analytical 
value for the chiral limit (as determined in ChPT) as the K ^ formula given in Eqs IMl 
The K TT result does not have a divergent part and the smallness of this amplitude in 
the numerical data compared to the K amplitude reflects this fact. The chiral limit 
agrees within statistical errors for both methods, which is better than expected, given the 
systematic errors associated with the K —* amplitude. 

The NLO contribution to the matrix element K tt was computed in order to control 
the systematic error associated with the extrapolation to the chiral limit. The linear fit is 
exact to NLO since the log term vanishes, and as Fig. d illustrates, this fit is quite good. 
We therefore conclude that the error in our value of <^ in lattice units is dominated by 
the statistical uncertainty, which is around 10%. As mentioned previously, these results are 
in rough agreement with the large estimate of Ref. Q], and so, as these authors pointed 
out, the term proportional to a^^ will have a large numerical effect on the matrix elements 
of Qe- 

VIII. CONCLUSION AND OUTLOOK 

The application of the quenched approximation to weak decay matrix elements is made 
more challenging by the combination of 1) the perturbative calculations (in which QCD 
vacuum polarization effects are important) needed to connect the high energy scale typical 
of the W, Z and top masses with the low energy scale of the actual decay and 2) the non- 



32 



perturbative QCD calculations (in which the quenched approximation might be employed) 
needed to evaluate the relevant low energy matrix elements. We have demonstrated in some 
detail how the quenched approximation can be applied to the latter without altering the 
calculations which underlie the former. The resulting approach is very close to that employed 
in the two large-scale quenched calculations reported in Refs. 

We then reviewed the quenched chiral perturbation theory results for the strong penguin 



amplitudes relevant for e'/e as presented by Golterman and Pallante 
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They have 



shown that a new LEG, a^^, contributes to quenched amplitudes of the strong penguin 
operators (specifically, to Qe) and they have proposed a lattice method to obtain this new 
constant. 

In Ref. a^^ was calculated using the large Nc approximation and the value was 
found to have a large effect on the K —>■ tttt matrix element of Qq. Such an effect could 
significantly alter the value of e'/e reported in Ref. Therefore, we considered it essential 
that a^^ be computed directly on the lattice. The method proposed by [3| to use the matrix 
element of ^ to obtain a^^ from the lattice was implemented in this paper, but it was 
shown that this method suffers from ambiguities due to power divergent contributions when 
evaluated using the domain wall fermion formalism. We have shown that the amplitude, 
K ^ IT, provides an alternative method to obtain a^^ from the lattice but without power 
divergent contributions. We implemented this method and obtained a value of a^^ which 
is indeed large enough to have an important effect on the K ^ im matrix element for Qq, 
and, therefore, also on the quenched determination of e'/ e. 

We conclude that we cannot reliably construct the matrix element {71^11^ \Qq\K^) within 
the quenched approximation using quenched ChPT and our lattice data. This is due to the 
large value of , a low energy constant absent outside of the quenched approximation, 
which was implicitly assumed to be zero in the previous RBC and CP-PACS work. Such a 
large quenching artifact does not merely make the extraction of {7i^7i^\Qq\K^) practically 
difficult, but implies that the quenched approximation itself fails to accurately describe the 
full theory. In fact, the motivation for the definition of the quenched approximation outlined 
in Section II is seen to be invalid, given the large differences in the analytic structure of the 
full and quenched theories. This clearly reduces the physical relevance of such a calculation. 

These arguments and results demonstrate that the errors associated with the quenched 
approximation in the evaluation of the strong penguin operators are likely to be quite sig- 
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nificant. Definitive answers will clearly have to await dynamical simulations. The first steps 
in performing such full QCD simulations have been underway for quite some time with two 
dynamical flavors ^3]; the more realistic 2+1 flavor simulations have begun. 
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APPENDIX A 



The function, F, appearing in the one-loop contribution to K ^ is given by 



' J p^ -\-m\ + m2— p^ a/ A(1, m\/p^, m\/p^) ' 

(Al) 



with 



A(a;, y, z) = {x — y + + Axy. 



APPENDIX B: NORMALIZATION CONVENTIONS 

In this appendix we specify the sign and normalization conventions used in this paper. 
This includes both the fields and operators used in chiral perturbation theory (where we 
follow in most cases the conventions used in the earlier work of Golterman and Pallante) 
and the corresponding quantities defined on the quark level where we follow the conventions 
used in Ref. |2|. 

Equations (23) and (24) define the relative signs of the pseudo scalar fields 0, the pseudo- 
fermion fields x ^'^^ the ghost fields in the sense that a specific pair of SU (3|3)l x SU (3|3)r 
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matrices, {Ul, Uji) will transform these fields in a determined fashion: 

S S' = Ul^uI. (B1) 

Note this relation between the fields 0, x ^i-nd is still somewhat abstract because the 6x6 
matrices Ul and Ur contain both commuting and anti-commuting numbers. However, this 
description will be sufficient for our purposes because these same matrices transform the 
quark and pseudo-quark fields q and q: 

^ (J) ^-(:)-^(:)^ 

Here Pr/l = (1 =t 7^)/2 and as in the text, Q contains the three flavors of quarks q and 
the three flavors of pseudo quarks q and belongs to the Cartesian product representation 
SUi3\3)L X SUi3\3)R. 

The absolute normalization and sign for the S field is determined by the lowest order 
effective chiral Lagrangian for quenched QCD: 

£g).^ = :^str{9^S^9^S} + strjMtS + S+M} (B3) 

(written in Minkowski-space) once we specify that the parameter Bq is real and positive 
with a magnitude chosen so that the 6x6 matrix M is the quark mass matrix appearing in 
the fundamental QCD Lagrangian written below. 

The Dirac fields q and q are conventional, with the Minkowski space Lagrangian given in 
terms of Q by 

C = QiYD^ Q - Q{M^Pl + MPr) Q . (B4) 

The connection between quark fields and the corresponding quantities in the effective chiral 
Lagrangian is given by the equation: 

{QlUQr\ = ^^,,. (B5) 

This equation is a generalization of the usual relation between S and the quark fields to 
include the new variables which appear in the quenched case. It is uniquely determined 
by the combined requirements of flavor covariance (the left and right-hand sides transform 
identically under the SU{3\3)l x SU{3\3)r flavor transformations of Eqs. ()B1|) and ()B2|) ) 
and consistency with the equation for the chiral condensate: 

{uRu,) = ^lnZ[M,Mt] = -^(Sn), (B6) 
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where we examine the condensate associated with the up quark, the first component of Q, 
and treat the matrix elements Mi i and Mj*^ as independent variables. Note, the flavor 
covariance of Eq. ()B5|1 is lost and the resulting equation invalid if the order of the factors 
{Qnjj and (Ql)* is reversed, since these 6-component fields are a mixture of commuting and 
anti-commuting quantities. 

Following standard conventions, we identify the field in Eq. (24) with meson fields 
according to: 



77- -'K^/V2 + r]/VQ K° 



(B7) 



y K- -2r//V6 J 

where the specific fields appearing in Eq. ()B7|) above destroy the corresponding mesons. 
Similarly the ghost field x of Eq. 24 can be written in an identical fashion if a tilde is added 
to each of the meson fields, e.g. vr^ is replaced by tt"^. With these conventions we can then 
examine specific components of Eq. ()B5|) to lowest order in chiral perturbation theory and 
obtain the useful relations: 

s^^'d = zBofK' (j = 2,^ = 3) , (B8) 
s-f^d = -iBofK° (j = 5,z = 3). (B9) 

This same identification can be made by an appropriate generalization of the results in 



Appendix A of Ref. 3| to the quenched case. 

Remaining consistent with Eq. ()B5|) . we can easily write down expressions relating quark- 
level operator states and chiral-level operator states. For a given meson operator in terms 
of the underlying quark fields, we have 

d^5s\0) = tBofK'lO) = zSo/|X°) , (BIO) 
dj5s\0) = -iBof{Ky\0) = -tBof\k') , (BU) 

where Eq. ()B10|) is the standard relation between states in terms of the quark-level and chiral 
operators. Eq. ()B11|) is found by a chiral rotation on Eq. ()B10|) to give the corresponding 
relation between these fermionic meson states. To be concrete, we also include here the rules 
needed for evaluating matrix elements of operators with these x fields. These are given by 
chiral rotations on the appropriate K° field rules, and found to be 

{0\K%x)\K%k)) = e"'^ (B12) 
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and 

(^°(k')|^°(k)) = 2EK{2Ttf6%k' - k) . (B13) 

Unlike the standard meson fields like K^, the ordering in these expressions is rather impor- 
tant, since K^{K^y = —{K^YK^. Additionally, as stated earlier in the text, the quark and 
ghost-quark propagators are equal flavor by flavor (when properly ordered), such that 

{0\d{x)d{y)\0) = {0\d{x)d{y)\0) = SF{x,y) , (B14) 

where {'y^D^ + m)SF{x,y) = 5^'^\x — y), for the case of continuum, Euclidean fermions. 

Next we apply this same approach to relate the various components of the four-quark 
operator related to Qq'^^ as they appear in chiral perturbation theory and at the quark 
level: 

Of^' = tr,, {(QL),"(7^)*(QL)n tr,, {str [Q'^iYYQn'^ } = «f (sASt) (B15) 



where the diagonal matrix A is defined in Eq. (21), a and b are color indices, tro is a trace 
over the (implied) Dirac indices, and we are exploiting flavor covariance and the conventions 



of Ref. 



(irp 

m 



We can use this expression to evaluate two cases of importance by using 
= 5a5j2 and (Ae)^^ = 5i^5j^. In other words, multiplying Eq. ()B15|) by either or \q 
and taking the supertrace over the 5*^7(313) indices, we get 

Qe^^^ = str [AeO'^A^T + h.c. = a^'sti{\Y^AT)) + h.c. , (BIG) 
X^qQns'] ^ ^ ^ ^ ^str(A6SASt) + h.c. , (B17) 



= str 



which are precisely the expressions we see in Eqs. (14) and (19) for Q^^^, and Eqs. (42) 
and (43) for Q^^^. 

APPENDIX C 

The NLO corrections to K ^ n are calculated from one-loop insertions of the operator 
Eq. Hni as well as local operators that begin at NLO: 



37 



(5f ^ = tr[A6L^]tr[E^ZEL'^] 
= tT[Xe{J:^A^,L^}] 
ar = tr[A6{StZS,S}] 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 



= tr[A6SUS]tr[5]. (C6) 

Each operator has associated with it an a priori unknown LEG which we call cf^, with a 
scale dependence similar to that of Eq | 



~r,NS/ \ ~r,NS/ \ , "^"g '/* i pl / 

Ci {f^2)=c^ (/ii) + — In— . (C7) 



where the authors 



The coefficient of the scale dependence, rji , can be determined as in 
apply background field and heat bath methods. The values for rji are 0, —2, —Nf/2, Nf/2, 
and 1 for i = 1 — 6 respectively. In this case Nf is the number of sea quarks, which for 
the quenched case should naively be taken to zero. However one must take care in this 
case because this operator fEq I43|) has been introduced precisely to add contractions of the 
four quark operators that contribute to the fermion determinant which would otherwise be 
absent in the quenched approximation in which we are working. 

We can determine the scale dependence in the quenched theory by making use of the 
following trick. In the partially quenched case, A Apg = 2(1 — 3/A^/, 1 — 3/A^/, 1 — 
3/Nf,—3/Nf,...,—3/Nf)diag, where the ffist 3 valence entries are 1 — 3/Nf, and the next 
Nf + 3 ghost and sea entries are —3/Nf. Taking the Nf — » limit in this matrix is singular, 
but when these factors multiply the rji factors above, this limit yields the correct scale 
dependence for the amplitudes in the quenched theory. 

The scale dependence obtained in this way agrees with that of a direct one-loop calculation 
of K ^ TT, which yields precisely zero in the case of degenerate quarks. Thus, the one-loop 
chiral log vanishes, and the only NLO contribution in the degenerate case is proportional to 
m1,j {rrif + mres in fits to DW fermions) with an unknown coefficient. The diagrams needed 
are given in Fig. IHl The wave-function renormalization vanishes for tt and K for this matrix 
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TABLE I: The ratio ((0|Q^^^|K) - {0\Q'^^^\K))/{0\sj5d\K) for each often nondegenerate pairs 
of quark masses. The chiral limit of this ratio {md 0) can be used to obtain a^^ . 

rus md = 0.01 nid = 0.02 rud = 0.03 nid = 0.04 

0.02 -2.898(13) x lO^^ 
0.03 -2.918(13) X 10-2 -5.010(13) x 10"^ 
0.04 -2.932(13) x 10"^ -5.020(12) x 10"^ -7.088(13) x lO'^ 
0.05 -2.941(13) X 10^2 _5.026(12) x 10"^ -7.093(13) x 10^^ _9. 141(14) x 10"^ 

TABLE II: The matrix element (vf |Q^^^|-f^) as a function oi mj, the single quark mass that 
appears in the degenerate kaon and pion states. The chiral limit of this data can be used to obtain 
using Eq. IHSl 

0.01 1.702(69) X 10-2 
0.02 2.293(68) x 10"^ 
0.03 2.910(70) X 10-2 
0.04 3.545(72) x 10^2 
0.05 4.199(74) x 10^2 
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FIG. 1: A — > TT"*" diagram in which the vacuum polarization loop is entirely contained within 
a high-momentum subgraph. 
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FIG. 2: A — > tt"*" diagram in which the vacuum polarization subgraph must appear entirely 
within a low-momentum subgraph. 
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FIG. 3: A tt"*" diagram in which the vacuum polarization subgraph must be wholly inside 

or outside of any high-momentum subgraph. 
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FIG. 4: A —>■ tt"*" diagram in which the vacuum polarization subgraph can be partially contained 
within a high-momentum subgraph. 
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FIG. 5: An order ag diagram which generates a non- vanishing ghost amphtude in a quenched 
theory from operators with no direct ghost quark couphng. 
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FIG. 6: The ratio {{0\Q'^^^\K) - {0\Q^^^\K))/{0\sj5d\K) versus (circles) and + m^e, 
(squares) for the ten nondegenerate quark masses. The hues are hnear fits to the data of the form 
r/o + ??imrf and r/o + ??i(mrf + rures)- 
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FIG. 8: The matrix elements ((0|Q^^^|i^) - {0\Q^^^\K))/{0\s-/5d\K) as a function of ma (big 
circles) and md + rrires (squares) and f"^ /{2a^^'^^){T\\Q^^^\K) (small circles) as a function of 
equal to the degenerate masses of the quarks appearing in the K meson. For comparison, we plot 
the K ^ ^ matrix elements above. The pre-factor of {tt\Q^^^\K) re-scales the matrix element so 
that the chiral limits of the two methods are expected to agree according to ChPT. This graph 
shows that the chiral limits do agree, thus yielding similar values of a^^ . 
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FIG. 9: Diagrams needed to evaluate the NLO amplitude K ^ n. NLO corrections include tree- 
level diagrams with insertion of the NLO weak vertices (crossed circle), one-loop diagrams with 
insertions of the LO weak vertices (small filled circles) and the O(p^) strong vertices (big filled 
circle). The lines represent the propagators of mesons comprised of valence and ghost quarks. 



